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Abstract. - The momentum signatures in nonperturbative multiphoton pair production for general 
elliptic polarization electric fields are investigated by employing the real-time Dirac-Heisenberg- 
Wigner formalism. For a linearly polarized electric field we find that the positions of the nodes 
in momenta spectra of created pairs depend only on the electric field frequency. The polarization 
of external fields could not only change the node structures or even make the nodes disappear 
but also change the thresholds of pair production. The momentum signatures associated to the 
node positions in which the even-number-photon pair creation process is forbid could be used 
to distinguish the orbital angular momentum of created pairs on the momenta spectra. These 
distinguishable momentum signatures could be relevant for providing the output information of 
created particles and also the input information of ultrashort laser pulses. 


Introduction. — Electron-positron (EP) pair creation in strong electric fields Bi 
is a fundamentally important exploration to probe the quantum vacuum. As a famous 
prediction of quantum electrodynamics, experimental observation of this nonperturbative 
tunneling pair production is still absent due to the fact that present available laser fields 
are far below the Schwinger critical electric field Ecr = rnf fe ~ 1.32 x 10^®V/cm, where 
m is the electron rest mass and — e is the electron charge (we use h = c = 1). However, 
with the rapid development of laser technology, the stronger fields will be achieved in near 
future through some laser facilities [4]-[^. This improves greatly the hopes to realize the 
observation of EP pair production from vacuum |7-12 . For long time two different pair 
production mechanisms are identified by the nonperturbative pair creation process (7 ^ 1 ) 
and the perturbative multiphoton pair production process (7 ^ 1 ), where 7 = miojeE^ is 
the well known Keldysh adiabatic parameter 13 , w and Eg are the frequency and strength 
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of external electric fields, respectively. Fortunately, the perturbative multiphoton process 
has been experimentally observed at the Stanford Linear Accelerator Center (SLAC) [l4) . 

Compared to cases of 7 <C 1 and 7 1, obviously the case of 7 ~ 0{1) is seldom 

studied because of the difficulty to get a simple asymptotic formulae theoretically. However, 
more and more evidence indicates that the pair production in this intermediate regime is an 
important and interesting topic because it contains the signatures of multiphoton process 
and nonperturbative mechanism. Many novel phenomena are discovered and deepen the 
understanding of pair creation mechanisms 15 18 . For instance, the dynamically assisted 


Schwinger mechanism [16[|17| by the superimposition of a strong and slowly varying elec¬ 
tric field (nonperturbative mechanism) and a weak and rapidly changing one (multiphoton 
process) can dramatically enhance the pair creation rate. Recently the effective mass in 
the nonperturbative multiphoton pair creation regime 18 has been pointed out which will 


not only explain the peak positions of momentum distributions but also likely be observed 
directly in the laboratory. 

In this Letter, our studies are focused on the momentum signatures in the nonpertur¬ 
bative multiphoton pair creation for general elliptical polarized electric fields. By using 
the real-time Dirac-Heisenberg-Wigner (DHW) formalism 19-21 , we first revisit the mo¬ 


menta spectra of created pairs in a linearly polarized field and exhibit the node structures 
in the momenta spectra 22 . The interference carpets is found, which appears also in the 
atomic above-threshold ionization problem 23 . By studying the positions of nodes, we find 


that they depend only on the external field frequency rather than the Keldysh parameter. 
Further, when the elliptic polarized fields are considered, the node structures will change. 
Even the nodes will disappear for large values of polarization. Finally, we make the re¬ 
lation clearer about the node positions associated to the absence of even-number-photon 
pair creation. This study, on one hand, is valuable to deepen the understanding of both 
nonperturbative mechanism and perturbative multiphoton process in arbitrary polarization 
fields. On the other hand, it is also helpful to understand other relevant physical processes, 
such as cosmological pair production 24 , heavy ion collisions 


28,29 


25 , and ionization of atoms 


and molecules |26[|27| , especially the above-ionization with elliptically polarized laser pulses 


Theoretical formalism. — Here we consider the field in an antinode of the standing 
wave formed by two counterpropagating laser pulses with appropriate polarizations. Since 
the spatial scales of the EP pair production are smaller than the spatial focusing scales 
of the laser pulse, the spatial effects are not significant and the magnetic field is absent. 
Therefore, we have a spatially homogeneous and time-dependent electric field. For the 
general polarization the field can be written as 


E(t) = 


An 




exp ( - ^) [cos(a;t 


(j>),6sm{ujt + (/'), 0 ]'' 


( 1 ) 


where Eq is the maximal field strength, r defines the pulse duration, lo is the laser frequency, 
4> is the carrier phase, and — 1 < <5 < 1 represents the polarization of the electric field, the 
factor vTT^ is used to ensure the same laser intensity for different polarization fields. For 
convenience, we set r = 100 /m, (j) = 0 and <5 > 0 throughout this paper. 

Our numerical results are based on the DHW formalism which has been used to study 
the vacuum pair production for different electric fields, refer to the references 19 -21 . Par¬ 
ticularly, Ref. showed that the DHW formalism would be recovered to the quantum 
Vlasov equation (QVE) if the spatially homogeneous and time-dependent electric fields are 
linearly polarized. However one of the advantages of the DHW formalism is that it can also 
solve more complex electric fields such as Q. In order to precisely obtain the momentum 
distribution function of created EP pairs / for the spatially homogeneous electric field Q, 
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we adopt the trick used in 30 and the DHW formalism is reduced to 


/ = 1/2 (E^'J'Wg, 

Wg = "HgWg + 2(1 - /)0(Ei. 


( 2 ) 


Here the dot denotes a total time derivative, cei = (m/n(p), p/r2(p), 0, 0)^, n(p) = (m^ + 
is the total energy of electrons with the kinetic momentum p = q — eA(t), q is the 
canonical momentum, A(l) is the vector potential; wg represents an auxiliary 9-component 

^ is a 10 X 9 matrix, ^ = (0 Ig) is a 9 x 10 matrix, and 

/ —ep-E^/w^(p) —2px —2m \ 

Hg = - 2 px 0 0 

y 2 (m^-f p • p^)/m 0 ^ J 


vector; = 


—p^/m 0 


is a 9 X 9 matrix. Thus, we can get the one-particle momentum distribution function /(q, t) 
by solving Eq. ([^ with the initial conditions /(q, —oo) = wg(q, —oo) = 0. 

Results. — By solving Eq. ll). we plot the momenta spectra of created EP pairs in 
Fig. 0 for linearly polarized electric field (<5 = 0) with different Keldysh parameters, from 
(a) to (f), 7 = 2, 3, 4, 0.5, 0.75, and 1. One can see node structures in the rings of momenta 
spectra of created particles by using the more realistic electric field Q with a Gaussian 
envelope. From Fig. Sf) where 7 = 1 , we find that the nodes are located at the positions 
where the momentum is an even multiple of half the frequency for an even-number photon 
pair creation (cf. the 8 -photon ring with 10 nodes), however, the nodes are located at the 
positions where is an odd multiple of half the laser frequency for an odd-number photon 
process (cf. the 7-photon ring with 8 nodes). This result can be explained via the factor 
[1 -I- (—l)”''"^® cos (4ga;/w • arctany)] (where s is the spin of created particles, s = 0 for 
bosons and s = 1/2 for fermions, n is the absorbed photon number), in the approximate 
analytical solution of the pair production probability for a sinusoidal field with an infinite 
period 31 32 . Furthermore, from the factor, one can see that the node positions depend 


on the Keldysh parameter. 

To carefully analyze our results shown in Fig. we plot the momentum distribution 
function of n-photon pair creation fniqx, + 00 ) as a function of momentum q^ {qy > 0 , qz = 
0) in Fig. Obviously, (a) indicates that the nodes (the minimum values) of even-number 
photon pair creation are still located at the positions where the values of qx are integral 
multiples of the frequency for different w and the fixed Eq. Furthermore, in (b), we find 
that for a given frequency the node positions will not be changed by the electric field strength. 
These results imply that the positions of the nodes in the momenta spectra depend only on 
the field frequency rather than the values of 7. This still holds true for a sinusoidal field 
without Gaussian envelope. Gonsidering this result and the effective mass signatures 18 


in nonperturbative multiphoton pair production, the approximate analytical expression of 
particle momentum distributions 


31 can be modified as 33 


2u}^ 


/n(q, +00) ~ ^li^w(q)[l -h (-1)”+^® cos^*] 

TT 

X (i(2H,-nis (q) 


(3) 


where w{q) = exp{-TrE„/Eo ■ [3(7) -H Q&i( 7 )( 9 y + 9 ^)/^^] -f 62(7)5^/7712} with 62(7) = 

-76'i(7) and 61(7) = 5(7) -h ^ g '{' r )/ 2 ; the prime represents d/dy, and 5(7) = | fg(l - 

v? )l/2/[l + (^y) 

2 ]^/^dM; Q is a parameter used to indue the contributions of the high order 
terms of momentum; /3* = 27r5a;/a; and r 2 mis(q) = (q^ -|- is the effective energy 

of electrons. The semianalytical results (dashed blue lines) are shown in Fig. (c) and 
compared with the DHW results (solid red lines). We can see that these two results are 
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Fig. 1: (color online) Momenta spectra of created EP pairs in the {q^, qy) plane for linearly polarized 
electric field. The upper row for Eq — 0.1i?cr and the lower row for Eq = 0.4i?cr. The panels in 
each row, from left to right, are for lj = 0.2m, 0.3m, and 0.4m, respectively. 
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Fig. 2: (color online) Momentum distribution function of n-photon pair creation fn{qx,+oo) as a 
function of momentum q^ {qy >0, qz = 0). (a) for n = 8 of Fig. 10 of Fig. 0 b), and 14 

of Fig. [^a), from top to bottom, (b) corresponds to 8-photon process with the fixed frequency 
Lj = 0.4m and Eq/Ect = 0.4, 0.3, 0.2, and 0.1, from top to bottom, (c) shows the comparisons 
between the numerical results (solid red lines) and the approximate analytical ones (dashed blue 
lines) with the parameter Q = 0.88 for n = 6, 7, and 8 in Fig. [^c), from top to bottom. 
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agreement with each other very well 33 . In addition, our result is valuable to precisely 


measure the frequency of external fields with or without envelopes by analyzing the node 
positions in momentum spectra. For instance, the distance between any two adjacent nodes 
in a ring along the qx direction is exactly the field frequency. Our findings can also explain 
why the spacing between the neighbouring peaks in the center of the longitudinal momentum 
distribution f{qx^+oo) is equal to the laser frequency (cf. Fig. 2 in Ref. 34 ). When the 
electric field parameters in 34 are chosen as E = O.lFlcrjT = 100/m and a = ujt = 5, one 
can find that 7 ^ 1 , and the pair production is belong to the nonpertubative multiphoton 
regime. Therefore, the center of momenta spectra will be split according to Eq. ([^. A 
similar result can also be seen in Fig. 0 <i) or (e). 

For arbitrary polarized electric fields, the momenta spectra of created pairs in the {qx, qy) 
plane with = 0 are shown in Fig. It is found that the node structures in the particle 
momenta spectra gradually disappear as the polarization increases, and the rings become 
continuous and uniform for the circular polarization 5 = 1. This is because the rapidly 
rotating electric helds are more isotropic in the {x, y) plane. Additionally, we find that the 
number of created particles in the center of momenta spectra decreases as the polarization 
increases, and reaches a minimum at 5 = 1. Meanwhile it is found that the rings shrink 
with the polarization increasing. For example, the radius of the ring corresponding to 6 - 
photon pair production (the smallest ring) is about 0.628m for 5 = 0, while about 0.568m 
for 5 = 1. These indicate that a large value of polarization 5 increases the thresholds of 
multiphoton pair creation for a fixed field frequency. Altho ugh the electric field s we used can 
give the same expression of the effective mass, m* = mi/l + e^EQ/{2m'^uj'^), for different 
polarizations according to 18 , the thresholds still grow with the increase of polarizations. 


This is beyond the effective mass model in 18 


In addition, we also present the momenta spectra for different polarizations in the {qx, qz) 
plane with qy = 0 and in the {qy,qz) plane with qx = 0, respectively, see Fig. In the 
{qx, qz) plane, one can see that for a very small value of 5, the node structures are pronounced 
and their positions can be determined approximately by Eq. However, when the value of 
5 becomes large, the field Ey will be very strong, then the node structures will be disturbed 
by the combined effects of the field Ex and Ey, and finally disappear for a large value 
of polarization. In the {qy,qz) plane, it shows that a small value of 5 corresponds to a 
small value of Ey, thus the node structures are not obvious except near qy = 0. With 
the polarization increasing, the larger rings are present, and a channel along the z-axis is 
gradually opened up near qy = 0. 

To study the relation between the node structures in momenta spectra and the absence 
of even-number photon pair production, the distribution function f{+oo) as a function of 
the laser frequency is plotted in Fig. By analyzing Fig. Ba), we find that for a linearly 
polarized field 5 = 0, once the momentum component along the external field is zero, namely 
qx = 0, the even-number photon pair creation will not occur. A similar result can also be 
seen for 5^0, i.e., the even-number photon pair creation vanishes at qx = qy = 0 , see 
Fig. Bt>) and (c). The former result can be explained quantitatively employing Eq. (B. 
When qx = 0, the expression becomes / oc [1 -|- (—1)”+^] for EP pair production. This 
shows that the momentum distribution function has a minimum for even-number photon 
pair creation, i.e., the even-number photon process vanishes. Another result can also be 
interpreted qualitatively applying Eq. <§■ However, since there are two electric fields Ex 
and Ey in the {x, y) plane, the momentum qx in B should be reasonably replaced by 
{Qx Then one can find that once qx = qy = 0, the even-number photon process will 

not appear. The above results indicate that the absence of even-number photon pair creation 
at = 0 for 5 = 0 and = gj, = 0 for 5 7 ^ 0 is just one example of the node structures 
whose positions can be well explained by Eq. B- ^ more physical interpretation of the 
absence of even-number photon pair creation in the case of zero momentum is according to 
the C-parity selection rule 22 . Based on this explanation, we can deepen the understanding 


of the node structures shown in Eig. by assuming that the orbital angular momentum 
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Fig. 3: (color online) Momenta spectra of created EP pairs for different polarizations. Each column 
from top to bottom corresponds to 5 = 0.1, 0.5, 0.9, and 1, respectively. Each row from left to right 
shows the momentum spectrum in the [qx, qy) plane, in the {qx,qz) plane, and in the {qy, q^) plane, 
respectively. Other electric field parameters are chosen as Eq — 0.1y/2EcT, and to = 0.4m. 

(0AM) of created pairs I is related to Qx/uj as I ~ kqx/ui kqx/i^n/X) ^ kXqx {k is an 
integer and A is the reduced laser wave length). For instance, when I is odd, the C-parity 
of EP pairs is 1, then the even-number photon process having the C-parity (—i)^*"^" _ 
is possible, while the odd-number photon process having the C-parity (—1)°'^'^ = —1 is 
prohibited. If the assumption turns out to be true, one can obtain the 0AM for the first 
time by analysing the particle momentum spectra. Moreover, Fig. |^b) also presents that 
for the electric field with a large value of 6, the resonance pair production corresponding 
to a large photon number can not occur. Particularly, when i5 = 1, only one-photon pair 
creation is present. It seems that a rapidly rotating electric field will restrain the large 
photon-number pair production for zero momentum case. 
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Fig. 4: (color online) Momentum distribution function /(+oo) as a function of electric field fre¬ 
quency uj. (a) for linearly polarized electric field with q = 0, 52 , = 0.5m, and qy = 0.5m. Other 
unmentioned momenta are zero, (b) for different polarizations with zero momentum. From top to 
bottom, 5 = 0, 0.5, 0.9, and 1, respectively, (c) for 5 = 0, 0.5, and 1, with = 0.5m, qy = q^ = 0. 


Conclusions. — We have investigated momentum signatures in nonperturbative mul¬ 
tiphoton pair production for general elliptic polarized electric fields by using the real-time 
DHW formalism. It is found that the node positions in momenta spectra depend only on 
the field frequency for linear polarization and will be changed by other polarizations. Fur¬ 
thermore, the thresholds of multiphoton pair creation grow with the increase of polarization 
under the same laser intensity. The momentum signatures not only possibly provide us 
with the 0AM of created particles but also can present some signatures of ultrashort laser 
pulses, such as the frequency. We believe that these phenomena are expected to be observed 
more easily in atom ionization for a general elliptical polarized laser pulses, since the basic 
physical picture in the pair production is similar to that in atomic ionization problems but 
obviously the laser field is available for the latter. 
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